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We introduce lattice models with explicit N=2 supersymmetry. In these interacting models, the 
supersymmetry generators Q yield the Hamiltonian H = {Q + , Q - } on any graph. The degrees of 
freedom can be described as either fermions with hard cores, or as quantum dimers. The Hamiltonian 
of our simplest model contains a hopping term and a repulsive potential, as well as the hard-core 
repulsion. We discuss these models from a variety of perspectives: using a fundamental relation 
with conformal field theory, via the Bethe ansatz, and using cohomology methods. The simplest 
model provides a manifestly-supersymmetric lattice regulator for the supersymmetric point of the 
massless 1+1-dimensional Thirring (Luttinger) model. We discuss the ground-state structure of this 
same model on more complicated graphs, including a 2- leg ladder, and discuss some generalizations. 



Supersymmetry is an exceptionally powerful theoreti- 
cal tool. It often allows exact computations in field the- 
ory and string theory, even when interactions are strong. 
In this paper we introduce M = 2 supersymmetric lattice 
models describing interacting fermions and monomer- 
dimer systems. We show that the continuum limit of the 
simplest of these models, defined on a one-dimensional 
lattice, is a well-known 1 + 1-dimcnsional quantum field 
theory with J\f = 2 superconformal symmetry. 

Our definition of M—2 supersymmetry is that the 
Hamiltonian H is built from two nilpotcnt fcrmionic gen- 
erators denoted Q + and Q~ = (Q + Y (l). It is 



H = {Q+,Q-}. 



(1) 



The fact that Q + and Q~ commute with H follows from 
the nilpotency (Q + ) 2 = (Q~) 2 = 0. Our models also 
have a fermion-number symmetry generated by F with 
[F, Q ± ] = ±Q ± . We shall show how, in at least some 
cases, this lattice supersymmetry extends to a space- 
time super (conformal) symmetry in the field theory de- 
scribing the continuum limit. Lattice models with a 
symmetry involving fermionic generators, such as the t- 
J model at J=±2t, are often called "supersymmetric" 
in the condensed-matter literature, but do not have a 
Hamiltonian of the form eq. (|l|). 

All eigenvalues E of the Hamiltonian eq. ([j]) satisfy 
E > 0. All eigenstates form either singlet or doublet 
representations of the supersymmetry algebra. All states 
\g) with E = must be singlets: Q+\g) = Q~\g) = §. 
Conversely, all singlets must have E — 0. All the other 
eigenstates of H can be decomposed into doublets under 
the supersymmetry, and conversely any doublet represen- 
tation is an eigenstate. This is simple to prove: a dou- 
blet consists of two states \s),Q + \s), where Q~\s) = 0. 



It follows from the definition of H and the nilpotency 
of Q ± that both of these states are eigenstates of H 
with the same eigenvalue. All eigenstates can be decom- 
posed into doublets: the four-dimensional representation 
{\s'},Q-\s'},Q + \s'),Q + Q-\s')) is reducible. Let 

\s)^W)-^rQ + Q-W) 

where E s > is defined by H\s') = E s \s'}. Then 
Q-\s) = 0, and (|s),Q+|s» and (Q-\s'), Q + Q \s')) form 
two irreducible doublets. 

The models we introduce can be defined on any lattice 
(or actually, any graph) in any dimension. The simplest 
model involves a single species of fermion c,, placed at 
any site i of the lattice. The fermion obeys the usual anti- 
commutator {cj, Cj} = Sij, and the operator F — c\ c i 
counts the number of fermions. We impose the restriction 
that the fermions have hard cores, meaning that fermions 
are not allowed on neighboring sites. A hard-core fermion 
is created by c\V<i> , where the projection operator 'P<i> 
requires all sites neighboring i to be empty: 
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j next to i 



(2) 



On this space of states, the supersymmetry operators are 
defined by 



(3) 



It is easy to verify that (Q+) 2 = (Q~) 2 = 0. The Hamil- 
tonian is therefore 



i j next to i * 



(4) 
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The first term in the Hamiltonian allows fermions to hop 
to neighboring sites, with the projectors maintaining the 
hard-core repulsion. The second term favors having more 
fermions, as long as they are more than two sites from 
each other. Thus one can view it as a repulsive potential 
for fermions, in addition to the hard core. 

There are two key questions to try to answer. The first 
is: what properties can be computed exactly using the 
supersymmetry? We have already noted the positive en- 
ergy and the pairing in the excited-state spectrum, but 
these are just the simplest consequences of the supersym- 
metry. The second question is: what (if any) field theory 
describes the model in the continuum limit? 

To illustrate the power of supersymmetry, we find the 
ground states for a chain of six sites and periodic bound- 
ary conditions. First, we count all the states. There 
is one state |0) with /=0, and six states q|0) with /=1, 
while because of the hard cores, there are nine states with 
f=2, and two with /=3. The vacuum obeys Q~|0) = 

and Q+\0) = Ei=i4l°). so (l°)><3 + l°)) make U P a dou- 
blet. The remaining five states with f=l are all an- 
nihilated by Q~ , and Q + acts non-trivially on them. 
There are thus five doublets with (/, / + 1) = (1,2). 
The states with /=3 are both annihilated by Q + , and 
Q~ acts non-trivially on both, giving two doublets with 
(/, / + 1) = (2,3). This accounts for all the states in 
the theory, except for two states with f=2. These two 
cannot form a doublet, because they have same fermion 
number. They therefore must be singlets, so there are 
two E=0 ground states in this theory, both with f=2. 
With a little more work, one finds that they have eigen- 
values exp(±i7r/3) under translation by one site. 

A basic quantity in a supersymmetric theory is the 
Witten index [Q 



W = ti [(-l) F e ^ H ] 



(5) 



Because the two states in a doublet have the same energy, 
their contribution to W cancels, leaving the trace only 
over ground states and W independent of (5. W can thus 
be found by evaluating (|J) in the — > limit, where 
all states contribute with weight (— 1) F . For example, 
for the six-site chain discussed above, we confirm that 
W = 1 — 6 + 9 — 2 = 2. For our model on the cube, one 
finds 1^ = 1- 8 + 16-8 + 2 = 3. 

Computing W for the model eq. (Q) on a general 
graph poses a fascinating combinatorial problem. Co- 
homology theory is a powerful tool to compute the num- 
ber of ground states at any fermion number, and there- 
fore also W. The supersymmetry generator Q + satisfies 
(Q + ) 2 = 0, and the E = ground states are precisely 
the states \s) that satisfy Q + \s) = and that cannot be 
written in the form \s) = Q + \s'). Those states form what 
is called the cohomology of the operator Q + and for its 
computation a variety of techniques are available. These 
include the 'spectral sequence' technique, which can be 



applied in this context as follows: one splits the lattice 
in two sublattices, with corresponding fermion number 
operators F\ and F2, so that F — F% + F^. Q + can 
also be split as Q± + Q% , so that Qf increases Fi by 
one. The two operators and Q\ are nilpotent and 
anti-commute. The first step in the spectral sequence 
is to compute the cohomology of Q\ . Q% becomes an 
operator acting on this cohomology, and the second step 
is to compute the cohomology of Q\ on this subspace. 
Often the process terminates here, and the result is the 
cohomology of Q + . In general the procedure continues 
for a finite number of additional steps, see e.g. Q for 
details. A similar procedure exists for a decomposition 
of a lattice into n sublattices, and in particular when 
every sublattice consists of a single site. Applying this 
to the iV-site periodic chain, with F\ consisting of every 
third site, we find ground states solely at fermion num- 
ber / = int((A + l)/3). For a chain with N = 3p with 
p integer, we find two ground states (so W = 2(— l) p ), 
while for N = 3p ± 1, there is a single ground state (so 
W = (~l) p ). 

We now address our questions in the one-dimensional 
case, where the Hamiltonian eq. (||) is on an A-site chain: 



N 

h = Y, ( c *Vi + 4+i c ) Pi+2 + Pi-iPi+i] (6) 
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where = 1 — ejej is the projector on a single site. We 
take periodic boundary conditions, so indices are defined 
mod N. The translation operator T commutes with both 
H and F, and its eigenvalues t satisfying t N — 1 char- 
acterize the eigenstates of H. The Hamiltonian eq. (||) 
resembles a lattice version of the Thirring model, a 1+1- 
dimensional field theory with a four-fcrmion interaction 
term. Below we will make the connection precise. 

Before we apply the Bethe ansatz, we indicate how 
the space of states Hn oi the model can be obtained by 
applying a systematic 'fmitization' procedure fl to the 
chiral spectrum of a specific J\f — 2 superconformal field 
theory (SCFT) with central charge c = 1. This construc- 
tion follows two steps. In the first step, the full chiral 
Hilbert space of the SCFT is written in a 'quasi-particle' 
basis, with the fundamental quasi-particles forming a 
supersymmetry doublet with charges (1/3,-2/3). In 
the second step, the momenta of the individual quasi- 
particles are constrained to a maximum value in the or- 
der of N, corresponding to a discretization of the space 
direction of the SCFT with spacing of the order 1/N. 
This then leads to a truncated or 'fmitized' partition sum 
QN(q,w), which is closely related to the partition sum 
Z N (q,w) = Tr{q^ lo ^ T w^" 3 ^), which keeps track of 
the eigenvalues of the translation operator T and the 
fermion number F of the lattice model eq. (||). This en- 
tire construction, to be detailed elsewhere, respects the 
supersymmetry. This gives a rationale for the existence of 
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supersymmetry generators on the space Hn, and it pro- 
vides an independent way of determining the quantum 
numbers t and / of the supersymmetric ground states. 
We remark that an analogous construction based on the 
spinon basis of the simplest SU (2)-invariant CFT natu- 
rally leads to the space of states of the Heisenberg and 
Haldane-Shastry models for N spin- 1/2 degrees of free- 
dom. Clearly, many generalizations are possible. 

Having discussed the model eq. (^|) from different 
points of view, we now give some of the results of a Bcthc 
ansatz computation Q. This computation lets us iden- 
tify the continuum limit of the theory. An eigenstate of 
H with / fermions is of the form 

^)=^^ 1 ,i 2) ...i / )cl i c] 2 ...ct / |Q) (7) 

where we order 1 < i\ < ii — 1 < 13 — 2 . . . . Bethe's 
ansatz for the eigenstates of J? is Js) 

tp(ii,i 2 , ...if) = ^ Ap MpfVpa" 1 ■ ■ ■ Vpf 1 ■ ( 8 ) 
p 

for some numbers {fii, . . . , /i/} and Ap\ the sum is over 
permutations P of the set (1,2,..., /). By construction, 
the translation operator T has eigenvalue t — n£=i ^k 

The next step is to find highest-weight states under the 
symmetries of the model. In Bethe's case, the symmetry 
is the 0(3) of the Heisenberg spin chain; here it is the 
supersymmetry. Here, Q - ^ = requires 



The Bethe equations are / coupled polynomial equa- 
tions of order N. They cannot be solved in closed form, 
and to make further progress, one usually needs to take 
N large. In our case, however, the supersymmetry al- 
lows us to derive more results from the Bethe ansatz 
for finite N. Precisely, we define Wk in terms of fik as 
Wk = {pk - q)/{qfJ-k - 1), where q = exp(-i7r/3). Then 
Baxter's Q-function || Q(w) = Y[i=A w ~ w k) has zeroes 
at w — Wk- Defining 7Z(w) = Q(w)(l + w) N ~f , we find 
that for the Wk giving the ground state, 



K(q- A w) = tq~ n K(w)+t 



1 q N n{q 2 w). 



(11) 



We derive an explicit expression for lZ(w) in the sequel 
H , but from ( p"lj ) directly we can rederive / and t for the 
ground state(s). When N = 3p with p an integer, there 
are non-trivial solutions to ([ll]) only when / = N/3 and 
t = (— 1) exp(±j7r/3). For N ^ 3p, one has only a single 
solution with / = int((iV + l)/3), and t = (-l)^ -1 . 

We can see heuristically why the ground state has 
/ = int((7V + l)/3). The potential term in (g) alone is 
minimized by the state with a fermion on every third site; 
adding any more fermions forces fermions to be two sites 
away and raises the energy. The hopping term alone also 
discourages fermions from being only two sites away, be- 
cause it has negative eigenvalues when fermions can hop 
to an adjacent site and back again, and the hard cores 
prevent this if there is another fermion two sites away. 
The state with a fermion on every third site, 



i _1 (Mfc) 



N-S - 



n 



[i k [ij + 1- ^k 



[ik/ij 



(J-j 



(9) 



for all k = 1 . . . /. These are very similar to the Bcthc 
equations for the antiferromagnetic XXZ spin chain at 
A = ±1/2 §. The only difference in (g) is in the left- 
hand-side, which in the XXZ case reads ([ik) N r, where 
t^I corresponds to twisted boundary conditions. 

Demanding that a state be an eigenstate of a continu- 
ous symmetry fixes all the free parameters in the Bethe 
ansatz. The miracle of the ansatz is that in the Heisen- 
berg and other integrable models, this state is also an 
eigenstate of the Hamiltonian. In our case, supersymme- 
try provides the miracle. Any Bethe ansatz state obeying 
Q~\s) = is either a singlet or part of a doublet, and so 
it must be an eigenstate of H . Its energy is 



E = N - 2/ + 

k=l 



1 

Mi H 



(10) 



The supersymmetry doublets appear naturally within the 
Bethe ansatz. If (/ii, . . . [if) satisfies the Bethe equations, 
then the set (1,/xi, . . .[if) also must satisfy the Bcthc 
equations, and moreover, both sets have the same energy 
E. It is straightforward to check that the states associ- 
ated with these two sets are related by = Q + (j) U) . 
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resembles a Neel state for a Heisenberg anti-ferromagnet. 
It is not an eigenstate of the Hamiltonian: the full ground 
state is disordered. However, like the Neel state, we ex- 
pect this state to be a part of the ground state. Our 
derivations of fas = int((7V+l)/3) confirm this intuition. 
This heuristic picture also gives the fermion numbers of 
the low-lying excited states. The excitations include de- 
fects (domain walls) in the Neel-like state, such as 



'00»0»00»00> 
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The fermion number of this configuration is just one 
higher than that of the Neel-like state, and it has three 
identical defects. Since defects can be moved arbitrarily 
far apart with no change in the potential, it is natural to 
treat each defect as a quasiparticle with charge 1/3. The 
existence of fractional charge in 1 + 1 dimensions is an 
old story; this was first discovered in field theory H. 

Finally, we give the field theory describing the contin- 
uum limit of (Q). When taking large, one can rewrite 
the Bethe equations in terms of densities of roots M, 
and then derive integral equations (known as thermody- 
namic Bethe ansatz equations) yielding the free energy. 
Our model has the same thermodynamic equations as 
the XXZ chain at A = 1/2, so the two models coincide 
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in the continuum limit. The continuum limit of the XXZ 
chain is described by the massless Thirring model 8j , or 
equivalently a free massless boson $ with action [B 



S = 



2.9 



Jdxdt [(9 t $) 2 - (9,$) 



The continuum limit of the A = 1/2 model has g = 2/3; 
this is the simplest field theory with N = (2, 2) super- 
conformal symmetry Q. The (2, 2) means that there are 
two left and two right-moving supersymmetries: in the 
continuum limit the fermion decomposes into left- and 
right-moving components over the Fermi sea. The boson 
also can be decoupled into left and right pieces, so that 
$ = $ L + $ R , while its dual $ = g($ L - The 
states of the field theory are given by the vertex oper- 
ators V m . n — exp(im$ + of conformal dimensions 
^l,r = {m ± gn) 2 /(Ag). The four components of the 
Dirac fermion in the Thirring model are V±i ,±1/2, while 
the supersymmetry generators are = V±i ± 3 / 2 and 
T3/2- In a finite size £, the lowest-energy 
state is in the Neveu-Schwarz sector, where the Thirring 
fermion has anti-periodic boundary conditions. This 
state has £Vs = — 7r/(6£). The lowest-energy states in 
the Ramond (periodic boundary conditions) sector are 
given by \±) R = Vo,±i/2|0},ys; both have energy zero. 
States in this conformal field theory can be built up by 
operating with the "spinons" ±1/2. 

Comparing this supcrconformal field theory with the 
lattice model on N — 3p sites, we identify our two E = 
ground states with the two Ramond vacua; all other 
states of the lattice model are in the Ramond sector as 
well. The U(l) quantum number m corresponds to / — 
N/3 (the fermion number relative to the ground state). 
The spinons here have charge ±1/3, so it is natural to 
identify these with the fractionally-charged excitations in 
the lattice model. 

One can define a quantum dimer model by placing the 
variables cj on the links instead of the sites of a lattice. 
The product in 'P<i> is then over links which meet the 
link i, and the projections are precisely such that over- 
lap ping dimers are avoided. On the N = 3p-sitc chain. 



symmetric models by including more projectors in (g), or 
by including several species of fermions. Consider a two- 
species model, with fermions labeled by + and — , and 
with the conditions that (i) a single site may not be oc- 
cupied by two particles and (ii) same-type particles may 
not occupy nearest neighbor sites. On a periodic chain 
with 4n sites, the Witten index of this model turns out 
to be W — 3, and we have strong indications that the 
continuum limit of this theory is the second model (at 
c = 3/2) of the series of Af=2 minimal superconformal 
field theories. A typical ground state pattern is 



and one recognizes the possibility of domain walls of 

charge ±1/2 (o -\ o) and neutral defects (+ o +). 

The +/— pattern indicates an Ising substructure in the 
model, in accord with the fact that the c = 3/2 Af = 2 
minimal model can be written in terms of a Majorana 
fermion and a free boson. 

We finally remark that in our models, the space of 
states is made up solely of fermions on which the super- 
symmetry acts non-linearly. (This does not preclude a 
linear realization on bosons and fermions in the contin- 
uum theory.) Such realizations of J\f = 2 lattice super- 
symmetry seem very different from other known realiza- 
tions, for example in lattice gauge theory |pT| . 
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